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We present an alternative local definition of the writhe of a self-avoiding closed loop 
which differs from the traditional non-local definition by an integer. When studying dy- 
namics this difference is immaterial. We employ a formula due to Aldinger, Klapper and 
Tabor for the change in writhe and propose a set of local, link preserving dynamics in an 
attempt to unravel some puzzles about actin. 
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Fuller's ubiquitously used relation between the link, twist, and writhe of a closed 
ribbon |IJ has played a role in a great deal of work on the equilibrium statistical mechanics 
of DNA @,[|,|4|] • While the traditional definition of writhe unambiguously depends on the 
shape of the DNA backbone, it suffers from being non-local - it depends on the entire 
conformation. This non-locality is merely a technical complication for equilibrium physics 
but is a severe limitation for the study of dynamics. In this letter we present an alternative, 
local expression for writhe which differs from the usual expression for writhe only by an 
integer. Since continuous dynamical evolution cannot smoothly change an integer, from 
the point of view of dynamics this new local expression is adequate. We will, in addition, 
propose a local, dynamical link conservation law for a closed curve. Though we will not 
derive this conservation law from first principles dynamics, using arguments based only on 
the locality of interactions along a twist-storing polymer such as DNA, we will argue that 
at some scale it is appropriate. 

The classic result relates the linking number Lk of the two sugar-phosphate backbones 
of DNA to two quantities, the twist Tw, which is a measure of the rate at which one 
backbone twists around the other, and the writhe Wr, which is a measure of how twisted 
in space the average backbone is. These three quantities are related by 

Lk=Tw+Wr. (1) 

The double helix may be described by the average backbone curve R(s) with unit tangent 
vector T(s) and by a unit vector U(s) perpendicular to T(s) which points from the R(s) 
to one of the two backbones. In this case twist and writhe are 

Tw= jds T(s) x U(s) -U(s) (2) 

and 

Wr=±l ds l ds > tR(«)-R(Q]-R(»)xRfr') 

4tt / J |R(a)-R(a')| 
where X(s) = 4— [|5) and s measures the arc- length along the curve. We will first suggest 
a different local expression for writhe Wr which we will then show differs from the usual 
definition of writhe Wr by an integer. While there is no new mathematics in this derivation 
it is useful to have a local expression which captures much of the information of writhe. A 
similar resulting expression was first used in the context of Fermi-Bose transmutation by 
Polyakov 0. 
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We consider a closed curve V parameterized by R(s) with unit tangent vector T(s) 
and choose ei(s) and e 2 (s) along the curve to form an orthonormal triad with T(s): 
{ei(s), e 2 (s), T(s)}. Parameterizing the "difference direction" U(s) in terms of ei and e2 
we have 

U(s) = cos6 v (s) ei(s) + sin0 a (s) e 2 (s). (4) 

The total twist as defined in (||) is the integral of a local "twist density" which we can 
compute in terms of 9 V and e^: 

Tw= J dsTw = jdsTx\J-t = jds {d 8 9u - e%d s e%} (5) 

The first term on the right hand side of (|5]) is necessarily 2n times an integer. What about 
the second term? If we imagine that e»(s) = ej(R(s)) are defined everywhere in space 
then we define: 

2tt Wr = I ds e^d s e% = I dRi efdng. (6) 



j Jr 

Since are only defined on the curve, derivatives in space <9ie 2 are not necessarily well- 
defined. However, since we are only interested in this derivative along the curve, we can 
make sense of @. This expression appears to be a local expression for something which 
we would like, considering (|5|), to interpret as writhe (times 2tt). At first glance adding 
the orthonormal triad appears similar to the "generalized Frenet-Serret" frame introduced 
by Goriely and Tabor 0. In the following, we will propose a dynamics that makes no 
reference to the additional vectors ei and e.2- 

Suppose the curve V is the boundary of a surface M (we can find such a surface for 
an unknotted closed curve). Since T e S 2 and 7Ti(5' 2 ) = we can find a continuous vector 
field T on M which is equal to unit tangent vector T on the boundary r ||. In figure 1 
we summarize the geometry and notation. Using Stoke's theorem we then have 



I dR % e^die^ = [ dS, [V x e^Ve°] l . 
Jr Jm 



(7) 



The integrand appearing in the integral over M has significance in both two and three 
dimensions. If ei and e 2 are local tangent vectors to a two-dimensional surface embedded 
in three dimensions, then the curl of the "connection" u = e" Ve^ is the Gaussian curvature 
of the surface while in three dimensions if ei and e 2 are perpendicular to, for instance, 
the nuclear spin in 3 He or the long nematic director in a biaxial nematic then the curl of 
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u> is the geometric realization of a topological fact: in a biaxial nematic +1 disclinations 
cannot escape into the third dimension as in uniaxial nematics [TO] . The curl of the uj is 



e ijk djU k = [V x e?Ve^ = ^e^e^T^d^ d k T\ (8) 



the famed Mermin-Ho relation pl |. In the following we will show that there is a relation 



between our proposed expression for writhe and the more traditional expression 



Nearly twenty years ago [JT2J] Fuller showed that the writhe Wr of a curve differed by 
an integer from the signed area swept out by the unit tangent vector T on the unit sphere. 
For completeness we review that result here. As before, let T(s) be the unit tangent vector 
to the curve R(s) and U(s) be the vector pointing along the ribbon with U(s)-T(s) = 0. 
T(s) traces out a curve on the unit sphere (the so-called "tangent indicatrix" ) . Let the 
tangent to this curve be t(s). Since U(s) is perpendicular to T(s), it lies in the tangent 
plane of the unit sphere at T(s). Applying the Gauss-Bonnet theorem [[13] to the region 
A enclosed by the curve on the unit sphere gives 

d 2 cr K + / ds k 9 = 2n (9) 
JdA 

where K is the Gaussian curvature and K g is the geodesic curvature, defined as the rate 
of change of t(s) in the tangent plane of the sphere. Locally parameterizing A by the 
orthonormal vectors ei(o"i, a 2) and 62(01, 02) (with e\(<j) x 62(a) = n(<7i), where n is the 
surface normal), we may write the unit vectors t and U as 



t(s) = cos 6>t(s) ei(s) + sin6>t(s) e 2 (s) 
U(s) = cos9u(s) ei(s) + sin du(s) e2(s) 



(10) 



where e^(s) = e^(cr(s)) and a(s) is the parameterization of the boundary in terms of the 
surface coordinates a. Since on dA the surface normal is T(s), the geodesic curvature 
is K g = Txt-t while the twist density Tw is defined as 2-kTw = TxU-U. Using the 
orthonormality of {ei, e2, T} we have 

Kg = d s 9 t (s) - e?<9 s e£, (11) 

and 

2nTw = d s 9u{s) - e%d s e%. (12) 



3 



Since the difference is 2-irTw — K g = d s (9 u — 6 t ) after one circuit about the curve the 
integral of the difference changes by an integer multiple of 2tc. Using (Q) and the fact that 
the Gaussian curvature of the unit sphere is K = 1, Fuller found that 

— + Tw = mod 1, (13) 

where A is the signed area of the region enclosed by T(s) on the unit sphere and Tw = 
J dsTw. Comparing this to Lk = Tw + Wr, we see that 

Wr= — modi, (14) 

which establishes Fuller's result. Upon using Stoke's theorem and the Mermin-Ho relation 
we have: 

2n Wr = j ds e1d s e^ = ^J dS t e ljk e af3l f a d 3 f^d k f^ (15) 

which is precisely the area swept out by the curve T on the unit sphere! Thus we see that 

Wr = Wr mod 1. (16) 

While the total amount of writhe including the integral part is important for calculating 
ground states of a particular twisted ribbon @,|],[| when considering changes in writhe 
the integer is less important. In particular to study dynamics a local "writhe density" is 
probably essential. 

To this end, a useful result [14|] for the change in writhe of a curve as a function of 
some deformation A 

d x Wr(X) = ±- £dsT(s,\)-[d x T(s,\)xd s T(s,\)}, (17) 

can be used. Note that this result follows from the preceding discussion: ( |T?D (multiplied 
by dX) is the differential in the area swept out by the tangents of two closed curves with 
tangent vectors T(s, 0) and T(s, dX). Since the twist is really the local torsional strain of 
the polymer, we denote it as 0(s, t) and then we have 

d t Lk= jds {d t n + d t T-[d s TxT}} + d t n (18) 

where n is the integer difference between Wr and Wr. Since continuous evolution cannot 
lead to discontinuous changes in the integer n and since link is conserved we have 

= jds {d t n + d t T-[d s TxT}}. (19) 



This conservation law need not be satisfied locally: the curve can twist in one place and 
writhe at some distant location to satisfy (|T9|). While mathematically consistent this is 
not a physically plausible effect. We expect that the linking number is conserved locally 
and only changes through the diffusion of some "link current" j. With this assumption we 
have the local conservation law 

d a j = d t n + d t T-[d s TxT], (20) 

which satisfies ([19]). Note that although we introduced the framing vectors they do not 
appear in the conservation law. The invariance of Wr with respect to the choice of only 
holds for a closed curve. Nonetheless the form of fl20|) suggests that for dynamics we may 
apply the conservation law to open curves . 

We consider now the relaxational Rouse dynamics of a closed, twist-storing polymer 
with a free energy 

F= \jds |A(d s T) 2 + Cfi 2 - A (T 2 -l)} (21) 

where A and C are the bend and twist elastic constants, respectively, and A is the Lagrange 
multiplier which enforces the constraint that T be a unit vector |[15|| . The dynamical 
equations for T and O are 

d t n = -rj-^- = -cr n n (22a) 

SF 

d t T = -r T — = AT T d 2 T + T T AT (226) 
oT 

= T 2 - 1, (22c) 
where F n and T T are dissipation constants. Using ( |2"2"| ) we may rewrite (|5l]) as 

d s j = -cr n n - at t t- [d s T x a s 2 T] . (23) 

We recognize the second term on the right-hand side of (|23| ) as k 2 (s)t(s) where k(s) is 
the curvature of the curve and r(s) is the torsion . Note that while usually the torsion is 
ill-defined for a curve with zero curvature, there is no problem here because of the explicit 
factor of k 2 . This formulation differs from the work in [|16|] by the introduction of a link 



density current which allows a local constraint. Thus we see that the "link current" obeys 

8 s j = -CT n O - AT t k 2 t. (24) 
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This form of the link current allows us to locally understand how link moves: link moves 
when there is a twist and when the curve is non-planar (i.e. t ^ 0). Assuming a Fick's 
law form for the current j = —Dd s Ck in terms of the link density Ck, we finally have 

Dd 2 s Ck = CT n fi + AT t k 2 t. (25) 

Note that Fick's law implies that if the link is uniformly spread along the polymer then 
the current is constant. 

As a simple application of ( p0|) we can consider a closed ribbon of radius r = k _1 lying 
in the xy-plane. If we consider relaxational dynamics around the curved ground state we 
see that in this case ([25]) will depend linearly on T through the torsion r (as opposed 
to quadratically if k were 0). We can consider the dynamics of the ribbon in extreme 
cases. The diffusion constant D can be large (infinite) or small (zero). Additionally the 
twist modulus C can be large (infinite) or small (zero). If the diffusion constant is infinite 
then fl20|) is not much of a constraint: any writhe or twist deformation can be quickly 
compensated for by small local link deformations. This is the limit in which dynamics is 
non-local and in which the only constraint is (|19|). 

The case of zero diffusion is more interesting and corresponds to "ultra-local" dynamics 
in which link is conserved point by point. In this case we have j = and so 

CT n O = -AT t k 2 t. (26) 

If the twist modulus is or small then it is easy to change O and this is hardly a constraint. 
However, if the twist modulus is large or infinite then (]2^) implies that if the curve starts 
planar then it remains planar since r = and k ^ 0. This is, of course, only true 
if the polymer has an everywhere non-zero curvature. We note that this observation 
may explain some puzzling experiments on actin, another double-helical polymer ||T7|,|18 



Though Lk = Tw+ Wr applies only to closed curves, we would expect the local dynamics 
of closed twist-storing polymers to be the same as that for similar open polymers. In those 
experiments on open actin strands there were two surprising observations. The first was 



that in [|17|| , at long wavelengths, the persistence length seems to grow from around 2nm 
to roughly 10/Um though it remained uniform at 2\im until some crossover wavenumber 
q c . The experiment was done by observing actin filaments oscillating in a sample volume. 
In order to distinguish in-plane motion from out of plane-motion, an actin filament which 
did not completely remain in focus for the duration of the sampling was not used in the 
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data analysis. It is rather surprising that a polymer would remain in a single plane for 
an extended period |19|. In the light of (|26| ) a possible explanation emerges: if an actin 
filament has zero average curvature, then, on average, (|20|) is not a linear constraint since 
k = 0. However if the longest wavelength undulation modes have not equilibrated then 
it is possible for (k) ^ 0. This could explain both effects. The long wavelength modes 
do not have good statistics which is why they appear to have an anomalous persistence 
length and, at the same time, the non-equilibration of these modes forces the torsionally 
stiff actin filament to remain in the plane for some time until link can diffuse in from the 
ends. 

As a mathematical aside for the cogniscenti we note that when continuing T into 
the region M there is an ambiguity since TC2(S 2 ) = Z. However, the integrand in (|T5| ) 
is the skyrmion density ||20| , which, when integrated over the entire region can differ by 
2-7T times an integer for interior continuations which are homotopically distinct. This is 
not a problem since there is already a mod 1 ambiguity between Wr and A/(2n). Finally, 
we might consider the case of a knotted curve. In this case the curve can be unknotted 
by adding to the curve small loops. It can be shown in this case that the addition of a 
small segment of curve will change Wr by an integer, which is therefore unimportant for 
dynamics. 

I would like to thank J. Kas, T.C. Lubensky, J. Marko, P. Nelson and D. Pettey for 
discussions. I especially thank T.R. Powers for many conversations and a critical reading 
of this manuscript. This work was supported through NSF Grant DMR94-23114. 
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Figure Captions 



Fig. 1. Summary of geometry and notation used. The curve of interest is V with capping 
surface M. The unit tangent vector T to the curve can be mapped to the surface of 
the unit sphere. The area swept out by the tangent vector on the tangent spherical 
map is A, while t is the tangent vector to the curve (the tangent indicatrix) lying 
on the surface of the sphere. 
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